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Abstract 

We study the long-wavelength effective description of two general classes of charged dilatonic 
(asymptotically flat) black p-branes including D/NS/M-branes in ten and eleven dimensional su¬ 
pergravity. In particular, we consider gravitational brane solutions in a hydrodynamic derivative 
expansion (to first order) for arbitrary dilaton coupling and for general brane and co-dimension 
and determine their effective electro-fluid-dynamic descriptions by exacting the characterizing 
transport coefficients. We also investigate the stability properties of the corresponding hydrody¬ 
namic systems by analyzing their response to small long-wavelength perturbations. For branes 
carrying unsmeared charge, we find that in a certain regime of parameter space there exists 
a branch of stable charged configurations. This is in accordance with the expectation that 
D/NS/M-branes have stable configurations, except for the D5, D6, and NS5. In contrast, we 
find that Maxwell charged brane configurations are Gregory-Laflamme unstable independently 
of the charge and, in particular, verify that smeared configurations of DO-branes are unstable. 
Finally, we provide a modification to the mapping presented in arxiv:1211.2815 and utilize it to 
provide a non-trivial cross-check on a certain subset of our transport coefficients with the results 
of arXiv:1110.2320. 
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Introduction 

Charged gravitational string and brane solutions play an important role in string theory. A large 
class of these objects has for a long time been known to be the supergravity incarnations of the 
fundamental string and Dp-branes [1]. Indeed, these spatially extended geometries source a set of 
gauge potentials and are solutions of (an appropriate) supergravity that are exactly interpreted 
as the classical coherent state of a (large) stack of fundamental strings/branes each charged under 
the corresponding potentials. This simple, yet profoundly deep, observation lies at the very heart 
of the AdS/CFT correspondence [2]. In general, the weakly coupled low-energy effective dynamics 
of the F-string and the Dp-branes is well-understood and conventionally captured by the open 
string sector, i.e., the Nambu-Goto and Dirac-Born-Infeld (DBI) action, respectively. However, 
at strong coupling, the dynamics is more appropriately captured by the closed string sector viz. 
supergravity. It is thus natural to ask how the various aspects of the low-energy DBI descriptions 
of the effective dynamics translate to the supergravity picture. The framework to answer these 
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questions is, in part, provided by relativistic hydrodynamics. Indeed, the existence of such a 
connection is not far fetched since hydrodynamics is the natural local generalization of the global 
thermodynamics of black branes. 

The effective description of gravity as a hydrodynamic theory is of course by no means new. 
The reader might be familiar with the membrane paradigm [3] where the effective black hole 
dynamics is formulated in terms of a non-relativistic fluid that lives on the (stretched) horizon 
of the black hole. Perhaps more familiar is the celebrated fluid/gravity correspondence [4] which 
emerges as the natural hydrodynamic limit of the AdS/CFT correspondence. In its original for¬ 
mulation, the long-wavelength effective description of A/” = 4 SYM is identified with a relativistic 
hydrodynamic theory. Here the transport coefficients (for a review on relativistic fluid mechan¬ 
ics, see [5]) are directly computed from gravity using linear response theory in [6] and from a 
direct perturbative gravitational computation in [4]. Moreover, various generalizations including 
charged cases, have been carried out in these schemes, see e.g. [7]. The latter approach is quite 
similar in spirit to the one considered in the present work and has previously been applied to the 
Schwarzschild black p-brane in [8] and for the Reissner-Nordstrom black brane in [9]. However, 
we emphasize that in this work we do not consider asymptotically (co-dimension 1) AdS branes 
but rather asymptotically flat conhgurations of general co-dimension. 

In more detail, in this paper we shall be interested in the effective behavior of a quite general 
class of black brane solutions captured by the action (1.1) which will be introduced in Sec. 1. 
Although this action is rather general, in this work, we will restrict ourselves to the cases where 
the black p-branes source a (p-|- l)-form or a 1-form gauge potential. In the following we will refer 
to these two types of brane charge as fundamental charge and Maxwell charge, respectively. In 
particular, the treatment of these two cases will allow us to capture the effective hydrodynamic 
descriptions of the NS and Dp-branes of type H supergravity along with the M2 and M5-branes 
of eleven dimensional supergravity. Moreover, our computation also captures the effective theory 
of (p -|- l)-dimensional smeared brane configurations of DO of type HA supergravity. However, 
instead of fixing the value of the dilaton coupling and spacetime dimension to the ones relevant 
for the specific supergravities, we will keep these parameters free. This allows us to extract the 
dependence of these parameters in the hydrodynamic transport coefficients, which turns out to be 
quite useful for examining the general features of the hydrodynamics. We note that the part of the 
computation pertaining to Maxwell charge also provides the non-trivial dilatonic generalization 
of the results for the Reissner-Nordstrom black brane worked out in [9]. 

Considerable research has gone into understanding the stability properties of black branes 
carrying various types of smeared charges (the literature is extensive, for a review see e.g. [10]). 
Understanding these properties is interesting from a pure general relativistic point of view viz. 
the Gregory-Laffamme (GL) instability [11], black hole phase transitions [12] etc., but also plays 
an important role for understanding aspects of the (un)stable vacua of string and M-theory. In 
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general, having access to an effective fluid dynamic description naturally allows one to address 
questions regarding the stability properties of the system in question. For example, to leading 
order, i.e., at the perfect fluid level (equivalently; at the thermodynamic level), an imaginary speed 
of sound signihes a fundamental instability in the system. In the context of brane physics, this was 
noted in [13] where an instability in the sound mode of the effective fluid of the Schwarzschild black 
brane exactly was identified with the GL instability of the brane. A GL instability is naturally 
characterized by a dispersion relation which describes the “dispersion” of the instability on the 
worldvolume. Although the exact form of this relation is only accessible numerically many of its 
features can be understood from the hydrodynamic description. In particular, turning the picture 
around, the lack of unstable hydrodynamic modes in the effective theory (to a given order in the 
derivative expansion) is equivalent to the lack of a GL instability, at least to that given order. 

One of the main results of this work is that (to hrst order in the fluid derivative expansion) the 
stability properties of the fundamentally charged dilatonic black brane remains solely determined 
by the speed of sound. In this way, the system will exhibit a stable branch of configurations for 
sufficiently small values of the dilaton coupling. This includes the Dp-branes, p < 5, of type II 
supergravity along with the M2 and M5 branes of M-theory. However, it does not include the 
D5, D6 and NS5 brane. We note that a similar conclusion was reached in [14], but we refine the 
analysis to Hrst order and use the correct value for the bulk viscosity. This result is furthermore 
in accordance with expectations from various numerical works [11, 15]. On the other hand, we 
find that the dilatonic Maxwell charged black branes cannot be made stable for any value of the 
dilaton coupling or the charge. In particular, this includes the smeared DO configurations of type 
IIA supergravity. Again, this is in accordance with general expectations [16]. We note that, as in 
Ref. [9], this is a genuine first-order derivative effect, meaning that the instability is not visible 
to leading order. 

Another interesting aspect of our computation relates to various proposed hydrodynamic 
bounds in gravity [6, 17, 18]. As mentioned above, we keep the dilaton coupling free in our com¬ 
putations. In addition to elucidating the (in)stability properties of the various branes this also 
provides us with an extra turnable parameter to examine the possible violation of the hydrody¬ 
namic bounds. As expected, we find that the shear viscosity-to-entropy ratio bound r//s > l/dvr 
is saturated for our entire class of gravitational solutions. In contrast we find that the (holo¬ 
graphic) bulk viscosity “bound” is violated for the entire class of fundamentally charged branes 
thus providing a simple example of its violation. 

Although the setup employed in this work is conceptually quite different from the hydro- 
dynamic limit of AdS/GFT (the fluid/gravity correspondence), some of the results of the two 
approaches can nevertheless be related. This is perhaps not too surprising since the latter should 
in principle be obtainable from the former in the near-horizon limit. The precise connection 
between the effective fluid dynamic theory of the asymptotically flat (non-dilatonic) D3-brane 
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and the fluid/gravity correspondence on AdSs was explored in Ref. [19]. Here the effective grav¬ 
itational dynamics of the D3-brane subject to Dirichlet boundary conditions at an appropriate 
cutoff surface was considered along the lines of ideas introduced in the paper(s) [20]. It was 
directly shown that the fluid/gravity correspondence and the hydrodynamic effective theory are 
obtained as the two (most interesting) extremes where the cutoff surface is taken to be located 
in the near-horizon throat region and at spatial infinity, respectively. Perhaps more surprisingly, 
the results of the fluid/gravity correspondence can also be related to the effective hydrodynamics 
of the (neutral) Schwarzschild black p-brane. This was shown in the paper(s) [21] where the 
authors managed to derive a mapping from a certain class of asymptotically AdS solutions to 
asymptotically flat solutions and vice versa. This class of solutions (on each side of the mapping, 
respectively) exactly includes the Ricci-flat and the AdS Schwarzschild black branes, respectively. 
More technically, the mapping is established by noting that the hydrodynamic sector on either 
side (of the mapping) is completely included in a reduced lower dimensional theory. Solutions 
to these two reduced theories can then be related by a simple analytical continuation in the 
dimension of the particular space on which the reduction is performed (a sphere and a torus, 
respectively). In this way, solutions on one side can be reduced, analytically continued and con¬ 
sistently uplifted to the other side and vice versa. This, conceptually quite simple, procedure 
neatly allowed [21] to derive the second-order fluid dynamic transport coefficients of the asymp¬ 
totically flat Schwarzschild black brane from the results of the (co-dimension one) AdS black 
brane in general dimensions [22] . In this paper we will present a modified version of the mapping 
which allows us to include Maxwell charge. 

The paper is organized as follows. In the Sec. 1 we introduce the leading order (seed) geome¬ 
tries around which we will consider hydrodynamic perturbations. Before presenting the first-order 
results of the fundamentally charged black brane and analyzing their stability, we briefly review 
the associated thermodynamics and explain how the perturbative procedure is carried out. In 
Sec. 2 we turn our attention to Maxwell charged black branes and work out the first-order trans¬ 
port coefficients and discuss the associated hydrodynamic stability properties. We conclude the 
paper by performing a non-trivial cross-check of our results with known results from AdS by 
introducing a modified version of the “AdS/Ricci-flat correspondence”. 


1 Hydrodynamics of black ]3-branes 


In this work we shall be interested in D dimensional p-brane solutions (i.e., solutions with horizon 
topology X of the following action^ [23] 


I = 


^^R*l-2d(/A*d(/>--^ 7'(q+2) A * J'(q+2)^ 


^We work in units where GNewton = t^. 


( 1 . 1 ) 
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Here ^{q+ 2 ) = dC'(q_|_i), q € I, are the dilaton weighted held strengths associated with the 

gauge potentials (^(q+i) and I denotes the collective set of gauge potentials in the theory. Notice 
that some of the forms can be of the same rank nonetheless they have different couplings aq to the 
dilaton which distinguishes them. The action (1.1) is quite general, however, for specihc choices 
of I and dilaton couplings, it captures the (bosonic part of the) supergravity descriptions of type 
IIA/B (in the Einstein frame) and M-theory relevant for description of the D/NS/M-branes and 
their intersections. Notice that the aforementioned supergravity actions also contain a topological 
term needed to preserve supersymmetry of the full theory, however, this term does not play a role 
for obtaining the flat p-brane solutions [24] . The action (1.1) corresponds to IIA (IIB) supergravity 
for D = 10, Xrr = {0,2} (2rr = {1,3}) and Xns = {1} [25] and eleven-dimensional supergravity 
for D = 11 and Xm = {2} [26]. In general, the D/NS/M-branes couple both electrically and 
magnetically to the above potentials. We unify the description in the standard way by writing 
the field strengths in the electric ansatz, where now the index q in the action (1.1) runs over the 
(allowed) spatial dimensions of the branes of the theory. Given a dilaton coupling a, it will be 
convenient to define a parameter N to further unify the description through the relation 

2_4 2{q + l){D-q-3) 

““AT D-2 • 


The real positive parameter N (usually an integer for string/M-theory corresponding to the num¬ 
ber of different types of branes in an intersection [27]) in general is preserved under dimensional 
reduction [28].^ Also note that the parameter N mods out the Z 2 reflection symmetry of the 
solution space a —)• —a. For all the fundamental D/NS/M-branes of type II string theory and 
M-theory iV = 1, and aop = (3 — p)/2 for the Dp-branes, while = 1 for the FI (api = —1) and 
NS5 brane (onss = !)• Also notice that the dilaton coupling (1.2) vanishes for D = 11, q = 2,5, 
consistent with the fact that M-theory contains no dilaton. 

In the following we will consider singly fundamentally charged p-branes, i.e., q = p- We 
therefore (consistently) truncate the action (1.1) to X = {p}. Moreover, we will keep the dilaton 
coupling a = Up free. One easily works out the equations of motion. 


□(/ 


a 

4(p + 2)! 

1 

2(p + l)! 


, d (e‘^^ * X) = 0 

{X-Xr^+(^2{dcpf- 


4(p + 2)! 


(1.3) 


where denotes the Einstein tensor and p, i' label the spacetime directions. As explained in the 
introduction, we are interested in p-brane solutions to the theory (1.1) characterized by horizon 
topology X where we break the di symmetries on (and thus implicitly breaking dt 

as well) perturbatively while maintaining the SO(n-|-2) symmetry on the transverse sphere (note 
that the total spacetime dimension D is related to the positive integer parameter n according to 


^Since we require the dilaton to be physical (i.e., we require > 0), the parameter N is bounded from above 


N < 


2{D-2) 


- (D-q-3)(q+l)- 
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n = D — p —3). These perturbations exactly capture the hydrodynamic sector of the black brane 
as we review below. ^ The family of leading order unperturbed (seed) p-branes solutions to the 
EOMs (1.3) was worked out in Ref. [28] and takes the form 

ds 2 = h-1^ [ - fuaUb dx“dx^ + Aab dx“dx^ + j , 

-, ^ (1-4) 

Here we have applied a general boost {uau'^ = —1) to the solution in the p + 1 worldvolume 
directions labeled by x“ = (t, x*). The tensor Aab = Vab + UaUb is the projector onto the directions 
parallel to the brane but orthogonal to while *1 = dx*^ A ... A dx^ denotes the induced volume 
form on the brane geometry. The two functions / and h are given by, 

/(r) = l-(^) , /i(r) = 1 + 70 . (1.5) 

Here tq parameterizes the horizon radius and 70 parameterizes the charge of the solution. 

According to fluid/gravity lore there is a one-to-one correspondence between the solutions to 
the EOMs (1.3) around the solution (1.4), and the relativistic Navier-Stokes equations, 

divT = 0 , d*j = 0. (1.6) 

Here T = is the effective stress tensor and j is the effective current which are matched order by 
order in a relativistic fluid derivative expansion. The effective stress tensor and current encompass 
the asymptotic data of the perturbed solution and the correspondence allows one to reconstruct 
the full gravitational solution (to any given order in the derivatives) from these asymptotic tensor 
structures. At lowest order, i.e., no derivatives and flat intrinsic geometry, the correspondence is 
trivial as it is non-dynamical. Indeed, computing the asymptotic stress tensor and current of the 
solution (1.4), one obtains the stress tensor and current [14] 

Tab = QUaUb + PAab , j = Q-k 1 ■ (1.7) 



Here g, P and Q denotes the energy (density), pressure and charge of the brane and can be 
computed from the Gibbs free energy G, which in turn is computed from the on-shell action and 
takes the form 


G = Vol(S’"+^) 



( 1 . 8 ) 


The temperature T and the chemical potential 4* are the intensive thermodynamic variables which 
do not depend on Newton’s constant and are easily written in terms of the parameters xq and 70 
using standard methods, 


r = 


n 




4 > = 


' iV7o 

1 + 70 


(1.9) 


^Breaking the SO(n + 2) would roughly correspond to elastic perturbations, see [29] and related works. 
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Notice that the entropy s and charge Q are conjugate to T and $ and are computed from G in the 
usual way. The energy density g and pressure P are then derived using the Gibbs-Duhem relation 
w = q + P = Ts and the defining identity G = —P — ^Q. At lowest order and flat intrinsic 
geometry, the correspondence between fluid dynamics and gravity is therefore just a convenient 
repackaging of black hole thermodynamics in terms of a relativistic (perfect) fluid. However, 
note that if one abandons the requirement of flat intrinsic geometry, the statement becomes an 
equivalence between gravity and perfect fluid dynamics on a curved p-submanifold (known as the 
blackfold approach, see [30] for a review), which is a non-trivial statement. In this work we keep 
the intrinsic geometry flat, or equivalently, we do not perturb the transverse sphere 

1.1 The perturbative expansion 

In order to carry out the perturbative procedure we proceed in the standard way. Here we will 
give a brief summary of the computation and refer to Appendix A for many of the details (also 
see the papers [4, 8, 9]). In order to ensure that the perturbative problem is well-posed, we need 
to cast the metric (1.4) into Eddington-Finkelstein (EF) form (i.e., a coordinate transformation 
x°‘ —>■ (T°‘{r) tailored so that |dr| = 0). In these coordinates the metric (1.4) takes the form 

ds'^ = h~P^^ —fuaUb + Aabjda^'da^ — 2h^Uada^dr + h^r'^ . ( 1 - 10 ) 

Transforming the coordinates of course also induces a transformation of the gauge held, however, 
one can easily show that in the case of Schwarzschild —>■ EF coordinates the transformation of the 
gauge held can be undone by a suitable gauge transformation. We therefore keep the form of the 
gauge held (1.4) consistent with our gauge choice introduced below. According to the ideas of 
huid/gravity outlined above, the solution (1.10) (along with the associated matter helds) is part 
of a larger class of solutions ds'j for which the parameters tt“, rg and 70 (collectively denoted by 
^ = (rt“,ro,7o)) are worldvolume huctuating functions. This solution, which reduces to (1.10) for 
constant rt“, rg and 7 g, is in general unknown, but can be probed perturbatively in a derivative 
expansion. Therefore for long-wavelength perturbations, 

dsj = ds^ -|- dsg -|- 0((9^) , Aj = A Aq 0(9^), (fj = (f> (pg 0{d^) . (1.11) 

Here ds^ is the geometry (1.10) expanded to hrst order in worldvolume derivatives and ds| denotes 
the hrst-order correction coming from the full solution ds^ (and similarly for the matter helds). 
The main purpose of this section is to compute the hrst-order expansion of the full solution and 
extract the hrst-order effective currents. The EOMs exhibit a large gauge redundancy. In order 
to simplify the computations it is convenient to choose a (consistent) gauge where the transverse 
components of {Ag) are taken to zero and furthermore {gg)rr = 0, {gg)nn = 0. The latter gauge 
choice allows us to reduce over the transverse sphere effectively leaving us with a p -|- 2 

dimensional problem. Notice that although the transverse drops out of the problem, the 


7 



parameter n still plays an important role as it will enter the various coefficients in the resulting 
set of equations. Choosing the appropriate ansatz for the perturbations (collectively denoted by 
V’9 = {9di and plugging them into the EOMs (1.3) produces two sets of qualitatively 

different equations, schematically of the form 

Constraint: Cg + 0(5^) = 0 , Dynamical: i^q = sa(r) + 0{d‘^) . (1-12) 

Here Cg is an operator acting on the derivatives of the intrinsic fields ^ and and are 
two first-order linear differential operators acting on the perturbations ijjg as a function of the 
radial coordinate r (and only r) and finally sg(r) is a source term (which is a rational function 
in r whose coefficients depend on d^). We note that in order to obtain the dynamical equations 
one has to explicitly use the constraint equations. The constraint equations are independent of 
the radial coordinate and, as the name suggests, put constraints on (i.e., relations between) the 
derivatives of the intrinsic fields. As expected, the constraint equations are found to be equivalent 
to the conservation equations of the leading order perfect fluid stress tensor and current (1.7). 
Notice that the conservation of j takes the form daQ = 0, which just expresses that the charge Q 
is non-fluctuating along the worldvolume. Locally the inverses of the differential operators 
and Lr ' exist and the formal solution to the dynamical equations reads 

i/jg = + (42))"^ hi+h2 . (1.13) 

Here hi and /i2 denote functions in the kernel of l 4 and respectively. We note that the 
inverse operators of and in general involve various integrations which when evaluated 
on SQ are quite complicated leading to various types of hypergeometric functions. The local 
decomposition (1.13) is of course not unique, however, if we require the inverse of l 4 to exist 
globally (in particular at r = ro) or equivalently require horizon regularity, the homogeneous 
solution hi is forced to vanish leaving only /12, which is then fixed by the boundary conditions 
and choice of fluid frame. Indeed, a subset of the freedom in the homogeneous solution /i2 
comes directly from the seed solution (1.10) and can be generated by 0{d) shifts vq ^ vq + 5rQ, 
7o 7o + <^70 and ^ -|- 6u°‘ (along with gauge transformations of Aj). Such shifts of course 

map to new (regular) asymptotically flat solutions which differ from (1.10) by 0(d). From a 
relativistic fluid point of view, this freedom is expected and corresponds to 0[d) redehnitions of 
the temperature, chemical potential and fluid velocity. In order to fix the gauge (the fluid frame), 
we require that the first-order solution reduces to (1-10) when the fluctuations vanish, which in 
turn corresponds to choosing the Landau frame on the fluid side. The remaining freedom in 
^2 parameterizes non-renormalizable modes which we require to vanish by virtue of asymptotic 
flatness.'^ In this way the freedom in the homogeneous solution /12 can be completely removed 
and we obtain the full solution to first order in the derivative expansion. 

"^For n = 1 one needs to fix some additional gauge freedom in ggri in order to ensure asymptotically flatness. 
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1.2 Transport coefficients 

Having obtained the first-order regular asymtotically flat corrected solution and reexpressing it 
in Schwarzschild coordinates, it is straightforward to compute the induced effective stress tensor 
using familiar techniques and extract the first-order transport coefficients. By direct computation 
we obtain the following effective stress tensor 

Tab = QUaUb + PAab + H^^^ 0(5^) with = -2r]aab - C'&Aab ■ (1.14) 


Here aab and 'd are the (fluid) shear and expansion of the congruences u°‘, respectively, and the 
coefficients r] and C are the corresponding shear and bulk viscosity.^ They explicitly evaluate to 


A ^ = 2 + (n + l)(l + (C-2n)7o) \ 

47r ’ T] n + l + Cyo {n + l + Cyo)^ ) 


(1.15) 


where we have defined the constant C = 2—n{N—2). For fixed tq, the viscosities are parametrized 
by the charge parameter 70 and the parameter N. The neutral limit can be obtained indepen¬ 
dently by taking either of the parameters to zero. Indeed, for large values of the dilaton coupling 
the dynamics of the brane effectively reduces to that of the neutral black p-brane (we have checked 
that our results reduce to those obtained in the neutral limit [8] ). We note that the shear viscosity 
increases with the charge while bulk-to-shear viscosity ratio decreases as charge is added to the 
brane (for fixed rg). As expected the gravitational solution saturates the hydrodynamic bound 
v/s P l/dvr [17]. We also note that the charge current j = Q -kl does not receive any first-order 
contribution. This is of course just a reflection of current conservation in the effective theory. 
Finally, we have checked that our results agree with those obtained for the D3-brane in Ref. [19] 
(here p = 3, D = n + p + 3 = 10, a = 0) 


1.3 Dynamical stability 

With the fluid transport determined, we can now address the stability properties of fundamentally 
charged branes by analyzing the response of their corresponding effective fluids under small long- 
wavelength perturbations. The dynamics of a charged fluid is governed by the worldvolume 
conservation equations given by Eq. (1.6). Here, we are interested in the case where the stress 
tensor T takes the form of a general dissipative fluid to first order (1.14). For the sound mode(s) 
we write the dispersion relation, valid up to second order in the wave vector k = as 

uj{k) = Pcsk + iask'^ . (1-16) 

Here Cg is the speed of sound given by = ^P/dg, where the thermodynamical quantity kept 
fixed when taking the derivative depends on the specific type of charged fluid in question. The 
^Notice that for p = 1 the shear tensor vanishes. 

®One has to employ the limit i?c —> 1 and make use of the identifications rl = ro 7 o, = ro(l + 70 ) and 
70 = (1 — Se)/Ss. 

^Here the relativistic wave vector is = {cj, A:*}. 
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attenuation coefficient controls the dampening of the longitudinal sound mode and is deter¬ 
mined by the first-order transport coefficients (1.15). In order for the fluid to be dynamically 
stable we must require that the speed of sound squared and the attenuation coefficient are 
both positive. In addition to the sound mode, the fluid also exhibits a transverse shear mode, 

ijj{k) = —k^ , ( 1 - 17 ) 

w 


which is stable provided that rj > 0. 

Now, for fundamentally charged p-branes the charge is not allowed to redistribute itself on 
the worldvolume of the brane, since it is conserved along all directions as noted previously. This 
means that the dynamics of the effective fluid will be reminiscent of that of a neutral brane, since 
the charge only plays a role in the equation of state. In particular, the leading order stability is 
solely determined by the sound mode (1.16). The system will therefore be stable to linear order 
in k as long as the speed of sound squared. 


f^'] = 1 + (2 - njV)7o 

yOgJ n + l + Cjo 


(1.18) 


is positive. We note that in the limit of vanishing charge (70 —>■ 0) we recover = -i/(n+i) [8] 
signifying that the neutral brane exhibits a Gregory-Laflamme (GL) instability [11]. However, 
for finite 70 we find, as in Ref. [14], that there exists a threshold. 


70 


1 

nN-2 ’ 


(1.19) 


above which the black brane is stable under long-wavelength perturbations to linear order. For 
fixed temperature, the threshold value 70 is precisely where the brane conhguration reaches its 
maximal charge. As illustrated in Fig. 1, this point exactly corresponds to a transition point 
between an unstable branch and a stable branch of configurations with the neutral brane con¬ 
hguration located at the endpoint 70 ^ 0 of the unstable branch. For sufficiently large dilaton 
coupling all fundamentally charged p-branes are therefore unstable, since they approximates the 
neutral brane. On the other hand, the stable branch connects with the extremal limit and since 70 
measures the ratio between local electrostatic energy and thermal energy of the huid, the stable 
regime is where the electrostatic energy is dominant. Notice that the limit 70 ^ 00 corresponds 
to hat space. However, it is not entirely clear if any of the hydrodynamic interpretation survives 
in the strict limit. The threshold value does not exist in instances where nN < 2, since 70 is 
a non-negative parameter. In those cases, the speed of sound is imaginary for all values of 70 . 
Finally, for the branes in ten and eleven dimensional supergravity we have listed several values 
of interest in Table 1. In particular, we note that the threshold 70 does not exist for p = 5,6 
in D = 10. This is in agreement with the expectation that the supergravity descriptions of 
D/NS/M-branes are stable with the exception of the D5, D 6 and NS5 brane [11, 14, 15]. 

We can now proceed by rehning the analysis to quadratic order in k by including the second- 
order term in the sound mode. As explained, in the absence of charge diffusion, the attenuation 
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Brane a 


70 


s/ri 


3-p 1 (5 -p)7o- 1 72-17p+p^+8(8-p)7o+4(9-p)7^ 

2 5-p 8-p+(9-p)7o (8-p+(9-p)7o)^ 


Dp 

M2 

M5 


470-1 

7+870 

70-1 

4+570 


1 + 


1 + 


63 

(7+870)^ 


(4+570)^ 


Table 1: The expansion coefficients of the sound mode (1-16) for the p-branes of ten and eleven 
dimensional supergravity {N = 1). For a sufficiently large value of 70 the threshold value 70 is 
exceeded and the fundamentally charged branes are stable. We note that the threshold value in 
ten dimensional supergravity increases with the spatial dimension and that black NS/Dp-branes 
are always unstable for p > 5. It is also worth noting that the values for the D1 and D4 are 
equivalent to those of the M2 and M5, respectively, due to Type II A ++ M-theory relation. 
Finally, the values for the fundamental string and NS5 are equivalent to those of the DI and D5, 
respectively, as N is invariant under a —)■ —a. 


coefficient will take the exact same form as a neutral fluid, hence 



The effects of the fluid carrying charge therefore only appears through the explicit dependence on 
Q in the shear and bulk viscosities given by Eq. (1.15). When the threshold (1-19) exists, we find, 
as shown in Fig. 1, that the attenuation coefficient is positive for both branches of configurations. 
The stability is therefore fully determined by the linear order, i.e., by the speed of sound. This 
is in contrast to configurations with smeared charges where the attenuation coefficient plays an 
important role for the stability properties of the effective fluid as we shall see in Sec. 2.2. When the 
threshold exists, the fundamentally charged p-brane is therefore dynamically stable for sufficiently 
large charge parameter 70 at least to next-to-leading order. 

An occurrence of a dynamical GL-like instability is conjectured to be intercorrelated with 
a thermodynamical stability [15, 31] viz. the correlated stability conjecture. It is therefore 
interesting to relate the above dynamical analysis to the thermodynamic properties of the effective 
fluid. The condition for thermodynamic stability is computed in the canonical ensemble, since the 
charge is fixed in the system and thus only requires positivity of the specific heat cq^. Using the 
thermodynamic quantities (1.9) and the expression for the speed of sound (1.18) one can show 
that the point where the specific heat becomes positive overlaps precisely with the threshold 
value given by Eq. (1.19), i.e., the configuration with the maximally allowed charge (for a given 
temperature). Indeed, as pointed out in Ref. [14], there is a direct relation between the speed of 
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Figure 1; The qualitative behavior of the coefficients of the sound mode given by Eq. (1.16) as a 
function of the charge Q for a fixed temperature T. The quantities and correspond to the 
speed of sound squared and the attenuation coefficient normalized with respect to their neutral 
values, respectively. There is one unstable branch connected to the neutral limit (70 —>■ 0) plotted 
with dashed lines and one stable branch connected to the extremal limit where both coefficients 
are positive. The stable branch is reached exactly when the charge parameter 70 exceeds the 
threshold 70 given by Eq. (1.19). This point corresponds to the maximal charge of the brane 
configuration for a given fixed temperature. 


sound and the specific heat cq^ at fixed charge given by 




dg) 


Qi 




s 

CQp 


( 1 . 21 ) 


Eor this system, we therefore find that the dynamical stability is in agreement with the correlated 
stability conjecture. Although, our results are only valid to first order we expect the statement 
to hold to all orders. 


2 Hydrodynamics of dilatonic Maxwell charged branes 

In this section, we write down the first-order effective hydrodynamic expansion for the p-branes 
of the theory (1.1) with g = 0 (Einstein-Maxwell-Dilaton (EMD) theory). These results provide 
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the dilaton generalization of the results originally presented in Ref. [9] . In terms of the parameter 
N, the dilaton coupling a now takes the form, 


a 


2 


N \D-2) 


( 2 . 1 ) 


In particular, for the parameters relevant for type IIA supergravity, we have a = 3/2, and the 
action (1.1) is that appropriate for describing smeared configurations of DO-branes. The leading 
order solution is given by 


= h ^ —/ttaUfedx^dx^ + ^Aafedx“dx^ + / ^dr^ + 


A 

= — log /i 


^(1) = • 


( 2 . 2 ) 


The effective stress tensor and Maxwell current are computed from the asymptotics of the solution 
and take the perfect fluid form, 


Tab — Q UaUb T P/^ab ; ja — Q ■ 


(2.3) 


Here the energy density g, pressure P and charge density Q are obtained from the free energy 
in the usual way. Direct computation of the on-shell action reveals that the free energy G again 
takes the from (1.8). Moreover, the two parameters vq and 70 are related to the temperature 
T and chemical potential $ as in Eq. (1.9). However, the Gibbs-Duhem relation now takes the 
form w = g + P = Ps + ^Q- Here we have used a calligraphed Q to denote the monopole charge 
density in order to distinguish it from the fundamental dipole type charge Q considered in the 
previous section. 


2.1 Transport coefficients 

Carrying out the perturbative computation of the first-order corrected fields and the correspond¬ 
ing effective currents follows the procedure explained in Sec. 1 (see Appendix A for many of the 
details and e.g. Ref. [9]). The most important difference between the two computations consists 
of the existence of an additional SO(p) vector (dynamical) equation in the Maxwell case. This 
is in accordance with the fact that a worldvolume one-form potential contains p more degrees of 
freedom than a top-form potential. The constraint equations coming from Einstein’s equations 
take the same form as before but the constraint equation deriving from the Maxwell equation 
is shown to be equivalent to current conservation, daj'^ = 0. This allows for fluid dynamical 
fluctuations in the charge density Q, which in turn shows up as a new charge diffusion trans¬ 
port coefficient in the effective theory. The most general first-order derivative corrected effective 
current, consistent with the second law of thermodynamics, takes the form (here written in the 
Landau frame) 

f = Qu- + Tl^^ + 0{d^) with = . (2.4) 
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Here £> is the transport coefficient associated with diffusion (appropriately normalized). 

Without further ado, we now present our results for the transport coefficients. The shear and 
bulk viscosities are given by 


s C 2 2 / (n + l)iV 

’'-S' + + 


(2.5) 


here C is the constant introduced below Eq. (1.15). Again, for fixed tq, the viscosities are 
parametrized by 70 and N and they reduce to the values of the neutral limit if either of the 
parameters are taken to zero. As expected, the shear viscosity saturates the bound rj/s > l/dvr. 
The diffusion constant D, associated with the effective current (Eq. (2.4)), is determined to 

® ^ 47rro(l+7o) 

V nA^7oA/(l + 7o)^' 

As is straightforward to check, these results agree with those of [9] in the limit where the dilaton 
coupling goes to zero, or equivalently, N —)• ■ 


2.2 Dynamical stability 

We now analyze the dispersion relations of the fluid associated to the Maxwell charged branes. In 
contrast to the fundamentally charged branes, discussed previously, the Maxwell charged branes 
exhibit charge diffusion which significantly changes the effective dynamics. Most notably, since the 
charge density can redistribute over the worldvolume, this gives rise to an additional longitudinal 
mode with the following dispersion relation, 

uj{k) = ia'^k^ . (2.7) 


This mode is associated to the attenuation of the (long-wavelength) diffusion of charge and is a 
first-order derivative effect. The existence of this charge diffusion mode explicitly means that the 
stability of the system to second order in k is not only determined by the attenuation coefficient 
Os of the sound mode (1.16), but also by the attenuation coefficient 053 . In order for the system 
to be stable to second order one must therefore require both to be positive. As before this system 
also exhibits a transverse shear mode given by Eq. (1.17). However, since the shear viscosity is 
positive, this mode does not play any role for the stability properties of the fluid. 

To linear order in k, the stability of the Maxwell system is, similarly to the fundamentally 
charged brane, solely determined by the sound mode. More precisely, it is dictated by the speed 
of sound squared which for the Maxwell system is given by 

= (?L] = ! + (2 8 ) 

* VSel . (l + 7 oW)(n + l + C 7 „) ■ ' ■ ' 


In the limit of vanishing charge (70 —>• 0) we again recover the neutral value for the speed of 
sound. Eor finite 70 we find the threshold. 


70 


1 


N-2 


(2.9) 
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Brane a 70 was/n 

DO i -1 - 8 -p+(9-p)-^o (9-^’)(8-p)(7o + 1)^c4 1 

Table 2: Coefficients of the longitudinal sound and diffusion mode for the black DO-brane smeared 
in p > 1 directions. Note that the hrst-order coefficient (speed of sound squared) is negative 
whereas both of the second-order attenuation coefficients are positive. 


above which the Maxwell charged black brane is stable to leading order. Note that the threshold 
only exists for > 2, or equivalently, when the dilaton coupling a is sufficiently small.® The 
qualitative behavior of the speed of sound (for N > 2) is illustrated in Fig. 2. In particular, we find 
that there exist two branches of brane configurations for the same charge density; an unstable 
branch (70 < 70 ) connected to the neutral brane configuration and a stable branch (70 > 70 ) 
connected to the extremal brane configuration. However, in contrast to the fundamentally charged 
brane, the merging point (defined by 70 ) does not coincide with the maximal charge configuration. 
On the other hand, if the threshold does not exist, i.e., N < 2, the speed of sound is imaginary 
for all charge densities. The system is thus unstable independent of the charge. This instability 
is in accordance with the one observed for configurations of smeared DO-branes {N = 1) [16]. In 
Table 2, we list the values of interest connected to this special case. 

We now address the stability to quadratic order in k. We therefore consider the two atten¬ 
uation coefficients and Oj) associated to the sound and diffusion mode, respectively. For the 
sound mode we find, in contrast to Eq. (1.20), the following modification due to the presence of 
charge diffusion, 


a-s 


1 

w 



^ C ^ 2 

^+2 + r ^ 


CqJ 



( 2 . 10 ) 


where we have introduced the specific heat capacity Cq and the (inverse) isothermal permittivity 
c given in Eq. (2.12). Eor the diffusion mode (2.7), we find 


2 ^ 

cjw Cq 


( 2 . 11 ) 


Assuming that N > 2, the attenuation coefficients exhibit a hyperbolic divergence around 70 . 
We refer to Fig. 2, where we have plotted a generic case of the two longitudinal modes. This 
divergent behavior at 70 is very different compared to the continuous behavior observed for in 
Sec. 1.3 (see Fig. 1). Indeed, it seems that the linearized analysis breaks down at the threshold 
7 o. However, besides this new feature we still find that both the speed of sound squared and 
the attenuation coefficient a* are positive when 70 surpasses the threshold. The stability of the 
sound mode is therefore, similarly to the fundamentally charged branes, solely dictated by the 
speed of sound. In contrast to the attenuation coefficient of the sound mode, it is quite apparent 
^Explicitly the dilaton coupling squared has to be smaller than . 
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Figure 2: The qualitative behavior of the coefficients of the longitudinal modes given by Eq. (1.16) 
and Eq. (2.7) as a function of the charge density Q for hxed local temperature T. The quantities 
c^, Us and as are given by the Eqs. (2.8), (2.10) and (2.11) normalized with respect to their 
individual neutral values. Dashed lines are for 70 < 70 while filled lines are for 70 > 70 . We 
observe that no region exists where all three quantities are positive at the same charge density 
Q. Furthermore, the attenuation coefficients both show a hyperbolic behavior at the threshold 70 
indicated by a vertical dashed line. Note that due to the choice of normalization of the diffusion 
coefficient in Eq. (2.4), the attenuation coefficient a® takes a fictitious finite value in the neutral 
limit. 


from Eq. (2.11), that the hyperbolic feature of the diffusion attenuation coefficient is simply 
dictated by the inverse speed of sound. We therefore find that the diffusion mode has exactly the 
complementary behavior of the sound mode, that is, when the sound mode is stable the charge 
diffusion mode is unstable and vice versa. The Maxwell charged brane configurations therefore 
suffer a GL instability for all charge densities. 

Finally, we check whether the dynamical GL-like instability observed above is connected with 
the thermodynamic stability of the system as predicted by the correlated stability conjecture. The 
conditions for thermodynamic stability of the Maxwell black brane are computed in the grand 
canonical ensemble since charge is allowed to redistribute itself in the directions of the brane. 
This exactly leads to the requirement of positive specific heat capacity and positive (inverse) 
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isothermal permittivity [32, 33], 


Cq = 

c = 



/ n + 1 + (770 \ 

V(niV-2)70-1;'" ’ 

( (70 + 1)(1 - [nN - 2 ) 70 )) sT 


( 2 . 12 ) 


It is straightforward to see that these two conditions Cq > 0 and c > 0 are complementary and 
can never be satisfied. Indeed, the two quantities exchange signs at 70 = ^/(nN- 2 ). Comparing 
with the dynamical analysis above, we observe that it is not sufficient to consider the dispersion 
relations to leading order, i.e., only considering the speed of sound c^, since it predicts the system 
to be stable above the threshold (2.9). However, when we include first-order corrections (the 
attenuation terms and aj)), we exactly identify a similar complementary behavior between the 
sound mode and the diffusion mode. We therefore find that the dynamical and thermodynamical 
analysis predict a similar behavior for the system, but point out that while the threshold value 
for the above thermodynamic quantities exactly corresponds to the point where the configuration 
obtains maximal charge density, the threshold (2.9) for the dynamical stability corresponds to a 
smaller charge density as noted above (with the coincidentally exception of n = 1). A similar 
effect was also observed in [9]. This can seem puzzling at first, but we emphasize that it is not 
in contradiction with the correlated stability conjecture. 


2.3 A check: Mapping to AdS 


In the previous section we have analyzed the hydrodynamic behavior of dilatonic Maxwell charged 
brane solutions. In this section, we present a non-trivial check of our results with Ref. [34] using 
a modified version of the so-called “AdS/Ricci-flat correspondence” [ 21 ] which in its original 
form provides a map relating a certain class of asymptotically locally AdS solutions to Ricci-flat 
solutions. Starting with pure Einstein gravity, we modify the map by performing a Kaluza-Klein 
(KK) reduction. This gives us the possibility to connect with EMD theory, which is a subset of 
the class of theories (1.1). Erom this starting point, we continue as Ref. [21] by performing a 
diagonal reduction and connect the reduced theory (after an appropriate analytic continuation) 
to theories which admit asymptotically local AdS solutions. With this modification in hand, we 
can map the first-order corrected solutions of the previous section to their corresponding class of 
asymptotically locally AdS spacetimes obtained in Ref. [34]. 

We begin by performing a KK reduction of the Einstein-Hilbert action. This takes us to 
EMD theory and connects with the action (1.1) in {np + 3)-dimensions with q = Q and fixes 
the coupling constant a to 


with a 
ansatz 


2 2(a + l) 

=- 


a 


(2.13) 


n + p + 1. We then perform a diagonal reduction (over the sphere S’"'”*'^) with metric 


ds2 = 



(2.14) 
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where the [p + 2)-dimensional reduced metric, the Maxwell field, the dilaton field, and the scalar 
field X s-re independent of the (n + l)-directions of sphere. Note that this ansatz includes the 
solutions given by Eq. (2.2). The action of the lower dimensional theory takes the form 


I = Yol(S^^^) [ e^f7^*l + *n(n + l)- ^ -- 

4+2 V «(a + 1) 


d(j) A'kd4> 


+ 


^— — —^dxA+dx—-e F 2 A-k F 2 ^ 


(2.15) 


where TZ is the [p + 2)-dimensional Ricci scalar and F 2 is the field strength of the KK gauge 
field. Note that the second term comes from the integration over the sphere and that we have 
performed a suitable Weyl rescaling 


2a{a + 1 ) 


(2.16) 


To establish the map, we now consider Einstein gravity with a negative cosmological constant 
A = —d{d — l)/2 in {d + l)-dimensions,® 


/a = 


/ d+1 


~ 2A ) * 1 


(2.17) 


We compactify the theory on a /3 = (d — p — l)-dimensional torus using the following ansatz 

" My-Aa(x,r)dx“j + • (2.18) 


d4 = dsl+2ix,r) + e 


Here y is a distinguished direction along the torus and d£|_^ denotes the line element of the 
remaining part of the torus. Using the reduction ansatz (2.18), the reduced action is (dropping a 
total derivative) 

Xa = Vo1(T^) [ e^lTZkl+kdid-l)- ^ d^Akd^ 

Jp+2 \ — 1 ) 

f d-l\ 1 2 ^ ^ 

+ ( —^ j dx A + dx--e PF 2 A-kF 2 j , (2.19) 

with F 2 = dA. Suppose that we now have access to closed analytic expressions of solutions to the 
reduced action (2.19) for any positive integer d. It then makes sense to extend the domain via 
an analytic continuation [35] of the parameter d. In this way, we can view solutions as a function 
of d which now can take any real value. In particular, it makes sense to consider solutions for 
negative values of d. The same reasoning can be applied to the action (2.15) where the solutions 
are viewed as functions of n. Direct inspection of the two actions (2.19) and (2.15) shows that 
they are proportional under the analytic continuation d ++ —n, 


VoI{S^+^)Ia ++ Vol(T^)X . 


( 2 . 20 ) 


This implies that we can obtain solutions to the reduced theory (2.15) by a simple reflection of 
the dimensionality parameter starting from the reduced theory (2.19) and vice versa. In this way, 
^Working in units where the AdS radius is set to unity. 
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knowing solutions of the form (2.14) of the higher dimensional theory (for any ri € N) allows us 
to consistently uplift to solutions of the form (2.18) (for any d — 1 € N) using the connection 
between the two reduced theories (2.20) and vice versa. 

Now, with the modified version of the mapping in hand, we can relate brane solutions of EMD 
theory + hydro to those of AdS^+i compactified on [d — p — 1) flat directions.The effective 
hydrodynamics of branes in AdS^+i compactihed on a torus (as described above) was worked out 
in the Ref. [34] directly from the AdS results in general (d +1) dimensions [22]. We can thus map 
the solutions of the reduced AdS theory to solutions of EMD considered in this work and in turn 
compare the transport coefficients. This provides us with a non-trivial cross-check of our results 
for the particular value (2.13) of the dilaton coupling. Going through this exercise, we find a 
perfect agreement between the transport coefficients (Eqs. (2.5)-(2.6) and Eqs. (3.4.38)-(3.4.40) 
of [34], respectively.^^). 

Finally, we can in a similar way as in Sec. 2.2 check the dynamical features of the asymp¬ 
totically local AdS solutions explicitly by computing the dispersion relations of the associated 
fluid. In contrast to the asymptotically flat solutions one finds that all the coefficients of both 
the sound mode (1.16) and the diffusion mode (2.7) are continuous and positive for all values 
of the charge density. We therefore find that no threshold value exists for the particular dilaton 
coupling (2.13) and furthermore that the hyperbolic behavior has been resolved through the an¬ 
alytic continuation (compare e.g. with Fig. 2). The solutions are thus stable to first order which 
is perhaps not too surprising, since they originate from stable neutral AdS branes [36] through 
consistent reductions. In addition, one can easily check that they are also thermodynamically 
stable by computing the specific heat and isothermal permittivity (see Eq. (2.12)) and checking 
explicitly that they are positive for all values of the charge. 

Discussion 

We have performed the perturbative procedure which captures the hydrodynamic sector of fun¬ 
damentally charged (dilatonic) black branes and branes with Maxwell charge smeared over their 
worldvolume. The main result of this paper are the first-order transport coefficients that deter¬ 
mine the dissipative behavior of the effective fluids. They are listed in Eq. (1.15) and Eqs. (2.5)- 
(2.6). Furthermore, for each class of branes we have obtained the dispersion relations of the 
effective fluids for which we refer to Fig. 1 and 2. As explained in the main text, many of 
our results apply to brane descriptions arising in string and M-theory. We have highlighted the 
interesting values in Table 1 and 2. 

We have computed the shear and bulk viscosities and in particular found their dependence 

^'^We note that care must be taken for fixing n = 1 after the analytic continuation, since while this might still 
be some perturbed solution it will not correspond to the asymptotically flat perturbed solution that we have been 
considering. In a similar way one should take care from the reversed perspective if fixing d = 2. 

^^Note, that [34] uses a different normalization of the diffusion coefficient. 
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on the (local) charge and dilaton coupling. We find that both classes of black branes have a 
shear viscosity satisfying r//s = l/dvr thus saturating the universality KSS bound [17]. It is 
reasonable to expect that this result holds for any kind of smeared brane solution of the action 
(1.1). Even though there, at least to our knowledge, does not exist a universal bound on the bulk- 
to-shear viscosity ratio, it is nevertheless still interesting to test our result against the well-known 
(holographic) “bound” proposed by Buchel [18]. We find that the ratio for fundamentally charged 
branes always violates this bound for all non-zero values of the charge and only saturates it in 
the neutral limit. In contrast, we find that branes with smeared charge continue to satisfy the 
bound for small charge parameter, but break the bound for sufficiently large values of the charge 
parameter. Although violations of the Buchel bound already exist in the literature [37], our setup 
thus provides two additional and physically simple examples. In Ref. [34] a different inequality 
and potential bound was proposed for Maxwell charged AdS branes, but it is always violated 
for our asymptotically flat branes. This was also observed in [38]. Indeed, if a universal bound 
on the shear-to-bulk exists, it would need a rather a non-trivial modification, it seems. Finally, 
the presence of smeared charge on the worldvolume gives rise to an additional hydrodynamic 
diffusion mode. We have determined the value of the transport coefficient D associated with 
charge diffusion and explicitly found its dependence on the dilaton coupling. In fact, we note 
that even though it generalizes the result for zero dilaton coupling in [9], it takes exactly the 
same form in terms of the parameter N. We note that all the transport coefficients are positive 
as required by thermodynamic consistency (the second law of thermodynamics). 

We have analyzed the dynamical stability of both classes of charged black branes by consider¬ 
ing the response of their corresponding effective fluids to small long-wavelength perturbations. To 
leading order, both systems have a branch of configurations that suffers from a GL instability in 
the sound mode [11]. This branch connects naturally with the neutral brane configuration which 
is indeed known to be unstable [8]. However, for a certain regime of the solution parameters, 
both systems also have a branch of stable configurations for a given charge connecting with the 
extremal limit. When taking first-order effects into account, this seemingly similar behavior of the 
two systems breaks down. For the system with fundamental charge we find that the attenuation 
coefficient of the sound mode is positive for all values of the charge and thus the stability is fully 
determined by the leading order i.e. by the speed of sound. In contrast, an additional longitudinal 
diffusion mode plays a crucial role for the Maxwell system which in turn leads to a complemen¬ 
tary behavior between the stability of the sound and diffusion mode. The Maxwell black brane 
therefore suffers from a GL instability for all values of the charge. This complementary behavior 
is similar to the behavior of thermodynamic stability under change of charge and hence our result 
is in accordance with the correlated stability conjecture [31]. We note that this was also observed 
for the Reissner-Nordstrom black brane considered in [9] and that similar observations have been 
done for other smeared branes, see e.g. [32, 39] for the case of smeared Dp-branes. 
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In Sec. 2.3, we presented a cross-check of the transport coefficients associated with the ffuid 
dynamics of the Maxwell charged branes utilizing a modified version of the “AdS/Ricci fiat 
correspondence” [21]. Although, it would be interesting to obtain a deeper understanding of 
whether any physics can be attributed to the existence of such a mapping, we can at least here 
comment with a partial answer with respect to its possible generalizations. In our endeavors trying 
to relate generic theories with an action of the type (1.1) to theories with a cosmological constant, 
we found that mappings of this kind are indeed possible. Unfortunately, despite allowing for a 
substantial amount of freedom between the field contents of the theories one in general ends up 
with pathological theories containing a scalar field having a kinetic term with the “wrong” sign. 
This is mainly due to severe restrictions arising from the requirements of consistent dimensional 
reductions. However, it is worth mentioning an interesting application of the current map (as 
it is presented in Sec. 2.3), namely, that it provides us with a tool for working out the second- 
order corrected charged solutions by starting directly from the known second-order results in AdS 
[4, 22]. We leave this exercise for the future. Before ending this paragraph, it is also interesting to 
note that since part of our computation is connected to the fluid description of the near horizon 
throat geometry (AdS) of stacks of D/M-branes, there are two, seemingly unrelated, schemes 
relating the results of the fluid/gravity correspondence to those of the asymptotically flat branes 
(and vice versa). It is therefore tempting to think that the two approaches are related although 
the precise connection still remains unclear. Hopefully our results will help shed some light over 
these unresolved issues. In this regard it is also interesting to note that the hydrodynamic regimes 
of the spinning black D3-brane was recently considered in [40] extending the work of [19]. Here, 
the theory naturally includes a U(l) gauge field. In the light of our results for the Maxwell black 
brane it would be of considerable interest to review whether they can be related (for specific 
parameters) to at least a subsector of the hydrodynamic limit of the spinning D3. 

In this paper we have been concerned with the hydrodynamic sector of fundamentally charged 
black brane solutions and provided the dilaton generalization of the results obtained for the 
Reissner-Nordstrom black brane [9]. In order to put these results in a broader perspective it 
would be interesting to investigate the hydrodynamic limit for branes with general smeared charge 
{q < p) and possibly multi-charged bound state solutions. These results will be presented in [41]. 
Furthermore, since smeared DO-brane solutions are related to charged Dp-branes by T-duality, it 
would be interesting to work out the Buscher rules in a derivative expansion necessary to relate 
the derivative expansion of the solutions. Quite remarkably, one would then be able to utilize 
the map we have provided to obtain the second-order fluid transport for all Dp-branes by simply 
starting from pure asymptotically local AdS solutions. We intent to pursue these ideas in the 
future. 
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A Details on the perturbative computation 

For the convenience of the diligent reader, we here provide some of the details on our perturbative 
computation outlined in Sec. 1.1. When solving the perturbative equations we focus on the 
Maxwell charged case as it is in many regards the most intricate due to extra dynamical freedom 
in the vector sector. In subsection A.3 we provide the most relevant equations (and differences) 
pertaining to the fundamentally charged case. 


A.l Setting up the perturbative problem 


As explained in Sec. 1.1, in order for the pertubative problem to be well-posed, we need to cast 
the fields into (ingoing) Eddington-Finkelstein (EF) form (a:“, r, D) —)• (a®, r, D) with (t“ = (u, cr*), 
defined in the usual way, i.e., so that |drp = 0. It is not difficult to verify that the coordinates. 




, r*(r) = r + J dr ^ 




(A.l) 


will do the job for both brane solutions (1.4) and (2.2). Note that we have chosen the integration 
constant appearing in r* so that r* —)• r for large r. In this way the EF coordinates reduce to 
ordinary Schwarzschild light cone coordinates far from the horizon. Notice that it is possible to 
write down a closed form expression for r* in terms of the hypergeometric Appell function Fi, 


r*(r) = rFi 


IN 1 , . 

— ;--,l;l--;l-h,l-/ 
n 2 n 


r 



1 

n — 1 r” 



(A.2) 


where the last equality applies for large r and is valid up to O ) ■ With this definition 

of r*, we will limit our analysis to the case where n > 2 (the analysis for n = 1 needs some 
modifications, however, at the end of the day, the results for the transport coefficients can be 
obtained by setting n = 1 in the results derived below). 
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In EF coordinates (A.l), the metric (1.4) for the fundamentally charged brane takes the form 


ds^ = h p+'i+i “ f'^aUb + d(7“d(T^ — 2h ^ Ua diT^dr + . (A.3) 

Similarly, the Maxwell charged brane (2.2) takes the form 

ds2 = /i"(^^)^^(^-/7Xa'Ub + /i^A„b)dcjW-2/if'«ad(T“dr + /i^r2dOj„^i)^ . (A.4) 

Transforming to the coordinates (A.l) introduces non-zero radial components to the two gauge 
fields. This transformation can, however, be undone by a suitable gauge transformation and the 
EF form of the gauge fields and in this particular gauge read 

A = --^(^) \/A^ 7 o( 7 o + 1)' * 1 and A = V^loilo + l)'uada°- , (A.5) 

for the fundamentally and Maxwell charged branes, respectively. In particular, the two gauge 
fields only have components in the brane directions. Also note that the dilaton remains invariant 
under the coordinate transformation as it is independent of the brane directions. 

We now promote the parameters tq and 70 to slowly varying worldvolume helds and look 
for the corrections ds^, Aq, (j)Q, so that dsp Af, (j)f solve the full set of EOMs to hrst order 
in the derivatives (cf. Eq. (1.11)). In order to do this, we first need to expand the leading 
order seed solutions (A.3)-(A.5) to hrst order in the derivatives. Carrying out this expansion 
is of course straightforward, however, the resulting expressions are rather lengthy and not very 
illuminating and therefore we omit them here. Moreover, in the following we employ the inherent 
Lorentz symmetry of the background to work in the rest frame in the point V around which we 
consider hydrodynamic huctuations (in the ultra-local sense). We therefore take |-p = (1, 0). 

Moreover, we dehne ro(cj“)|p = ro and 7 o(o'“)|p = 70 . Similarly, all the (worldvolume) derivatives 
are understood to be evaluated at V. For example, (9aro(cr“)|'p = daV^ and so forth. 

The rest frame has a residual SO{p) invariance which we use to split the resulting equations up 
into different sectors characterized according to their transformation under SO(p). For the fun¬ 
damentally charged black brane the scalar sector contains five scalars, the vector sector contains 
two vectors and the tensor sector contains one tensor. We parameterize the SO{p) perturbations 
according to 

Scalar: (Aa)„q...ip = -^^ 70(1 + 70 )', {ga)vr = h , 

-j_ Nn _ Nn 

{gd)vv = h , Tr {g 3 )ij = h -+p+iTr/ij , (I)q = , 

Nn \[f 1 \ 

Vector: {gg)^i = h "+p+^ fvi , {ga)ri = h V2 , (A.6) 

_ Nn _ 

Tensor: {ga)ij = h "+p+i , 

where ( 5 ^)^ = {ga)ij - ^{Ti{ga)u)Sij, i.e., the traceless part of {ga)ij and = fij - ^{Trfki)6ij. 
Similarly, the scalar sector of the Maxwell charged system contains five scalars while the vector 
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and tensor sector contain three vectors and one tensor, respectively. Here we parameterize the 
SO{p) perturbations in the three sectors according to 

Scalar: = -\/iV7o(l + 7o)' > {9a)vr = h V 2 , 

{gd)vv = h , Tr {gQ)ij = /i"+p+iTr/y , 4>d = U , 

Vector: {Aq)^ = - 1 /^ 70(1 + 70 ) a* , {gd)vi = /i"+p+i fvi , (A.7) 

{gd)Ti = /l^("2 + n+p+l)y^. ^ 

N _ 

Tensor: = /i"+p+i f- - . 

As explained in the outline of Sec. 1.1, we will work in a gauge where all the rest of the components 
of ds| and Aq are taken to be zero. The consistency of this gauge choice can be checked a 
posteriori. We now proceed with solving the resulting EOMs to hrst order in the derivatives. The 
EOMs (from the action (1.1)) take the form, 


Rfiu ~ fj,4A7u4> ~ — £^u + Oid"^) — 0 , 

W (c“*.P7...p,) ^ + 0(fl") = 0 , (A.8) 

^ +) + 0(8^) = 0 . 

Here g = 0 for the Maxwell system and q = p for the fundamentally charged system and T 
denotes the dilaton-weighted field strength T = e“‘^dA. Moreover, we have defined 


5 - 1 

‘^'^"“ 2(9 + 1 )! 


T7 T7 PO“-Pq 

^fipQ...Pq^U 


g + 1 

{D-2){q + 2) 



(A.9) 


Notice that the right-hand sides of the EOMs (A. 8 ), £’^ 1 /, and are all 0{d). The 

parameterization of the anzatze (A. 6 ) and (A.7) are exactly chosen in such a way that the resulting 
equations only contain derivatives of fab, ao...g and and will thus be directly integrable as 
explained below Eq. (1.13). In the two subsequent sections we provide the most important 
details for solving the two systems. 


A.2 Solving the Maxwell system 

In this section we give some of the details for solving the three SO(p) sectors for the Maxwell 
charged system. 

Scalars of SO(p) 

The scalar sector consists of eight independent equations which correspond to the vanishing of 
the components: £yy,£rv,£rr,'~£^£ij,£nn,£{(j)),A4y and Mr- (cf. Eq. (A. 8 )). 
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Constraint eqnations: There are two constraint equations; = 0 and A4r = 0. The two 
equations are solved consistently by 


9„ro = - 


ro(l - {N - 2)70) 
n + 1 + (2 - n{N - 2))7o 


diU^ 


dvlo = - 


270(1 + 7o) 


n + 1 + (2 - n(iV - 2))7 o 




(A.IO) 


The first equation corresponds to conservation of energy while the second equation can be inter¬ 
preted as current conservation. These are equivalent to the scalar conservation equations given 
by (1.6) in the rest frame. Under the assumption that the fluid configuration satisfies the above 
constraints one is left with six dynamical equations with five unknowns. 


Dynamical equations: The coupled system constituted by the dynamical equations is quite 
intractable. One approach to obtaining the solution to the system is to decouple the trace function 
Tifij. Once Tr/jj is known, it turns out, as will be presented below, that all the other functions 
can be obtained while ensuring that they are regular on the horizon. It is possible to obtain a 
3rd order ODE for Trfij by decoupling it through a number of steps. However, first it is useful 
to note that the particular combination of TrSij and leads to the equation 

A [r"+V(r)r'(r)] = -(9,n*)r- (2(n + 1) + C^7o) h{r)^-^ , (A.ll) 

where we have defined the constant (7 = 2 — n{N — 2) and 

T(r) = . (A.12) 

As we shall see this equation is very reminiscent of the equations for the tensor perturbations for 
which we know the solution to be 

T(r) = - 2((9iu*) ^ (1-F7 o)^ log/(r)^ . (A.13) 

Here we have imposed horizon regularity, since Tr fij and fj, are individually regular on the 
horizon. Once fj, is known in terms of Trfij we can use Srr to eliminate f^^ and then take 
linear combinations of the remaining equations. The resulting combinations can then be used to 
eliminate f^^ and f"y such that one is left with two equations in terms of a„ and Trfij which can 
then be decoupled by standard means. The resulting equation is schematically of the form 

[Tr/,,]'" (r) + [Tr/i,]" (r) + [Trfij]' (r) = STrir) , (A.14) 

where Hi,H 2 and do not depend on the sources (world-volume derivatives) and the source 
term Stt only depends on the scalar diU^. The expressions for these functions are however very 
long and have therefore been omitted. After some work, one finds that the equation is solved by 

Trfijir) = + 7oc^^ G{r) - 2{diu')Trf^f (r) , (A.15) 

where the terms containing the two integration constants and correspond to the homoge¬ 
neous solution. The entire family of homogeneous solutions to equation (A.14) of course has an 
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additional one-parameter freedom which has been absorbed in the particular solution Tvf^j\r) 
and been used to ensure horizon regularity. With the introduction of we can safely take 
= 0. The function G is given by 


G(r) = -- 


N 


2 + 2 - 


Nn 

—^-TT ) 

n + p 1J r"- 


-1 


(A.16) 


u + p + 1 r"- 

and has an intricate relation to the gauge choice {gd)Qn = 0 as we shall see later. The particular 
solution which is regular on the horizon is given by 


= ^(1 + 7o)^a7oG(?') + ~ +T)) " ) (1 + /37oG(r)) 


with the coefficients 


a = 2p 


2 (n + 1 ) + C* 7 o 


and P = p 


re + 2 + C* 7 o 


(re + 1)2 + C'7o(2(re + 1) + (770)7 ' ^\n + l + C^Q^ 

With Tr/jj (and f^) given, the equation = 0 will provide the derivative of frv, 


(A.17) 


(A.18) 


fwir) = 


2(re + 1) + (7^70 h{r) 2 


\+-i V dr - 


JS 


-'‘“s 




(A.19) 

Since this equation is a 1st order ODE, the regularity of the horizon is ensured by Tr/jj. Note that 
it is possible to perform integration by parts and use that the derivative of r* takes a simpler form. 
One can thereafter obtain an analytical expression for the resulting integral. This expression is, 
however, quite cumbersome and we therefore only provide its large r asymptotics 


nk r 
I i 


frvir) ~ + {dy) ^ 


(A. 20 ) 


k=l 


Here the homogeneous solution takes the form 


++) = c™ + 2p(n + p+l) + (» + rt(2p + C)g,„ 


2 p 2 (re+p+ 1 ) + ( 2 p + (7)p^7o 


(A. 21 ) 


and the particular solution is given in terms of the coefficients UrJ and which depend on 
n,p,a, and 70 . These coefficients are in general very long and their expressions are therefore 
omitted. 

Using the expression for in terms of Tr/^, the Maxwell equation = 0 becomes a 2nd 
order ODE for the gauge field perturbation, 

.2 A 


dr 


1 


I 


—a'v{r) 


nr 


2 (re + 1 ) + ( 7 -i: 7 o 


dr 


^[Tr^]'(r) + 4a^^^4(r) 


(A. 22 ) 


^^Note that equation (A. 14) has been derived under the assumption that diU^ + 0. This especially means that 
when there are no sources the one-parameter freedom disappears in accordance with (A. 15). 
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This equation is solved by a double integration. The inner integral is manifestly regular at the 
horizon, one can therefore work directly with the asymptotic behavior of the right-hand side 
before performing the integrations. The large r behavior of the perturbation function is thus 
found to be 






(r) -h {diu") - 


n 


r.nk 


n — 1 


+ 




k=l 


where the first term constitute the homogeneous solution, 


(A.23) 




= _L c' 


(r) = d 


(2) J2)!l 


70% 


‘Ip -|- C 


2p 


2{n+p + l) + ( 2 p-k (7)^70, 


(A.24) 


and the particular solution is given in terms of the coefficients and depending on n,p, o, 
and 7 o. 

The last perturbation function can be obtained from Tr£^jj = 0 which provides a 1st order 
ODE for the perturbation. Horizon regularity is therefore ensured by the horizon regularity of 
Tr/jj. Using the expression for in terms of Tr/jj the equation is schematically of the form 


f'vvi'^) = Gi[^^fij{r)] + G2[av{r)\+G'i[f^{r)] + Sii{r) , (A.25) 


where Gi,G 2 ,G 3 are differential operators and the source Su depends on Again, the full 

expressions have been omitted and we only provide the large r behavior, 


fvvir) « + (diu^) ^ 


„nk 

_0_ 

rUk 


k=l 


aivr + 


(A.26) 


with the homogeneous part given by 


fiv\r) =c^vv + + 7o)(cl^^ - 


( 1 ) (n-kp-%)(l-F 7 o)a^ 

2p 


^ ( {n + p+l + {2p + C)'yo)h{r) - p'^oNfjr) \ ^( 2 ) 

V p (2{n + p + 1) + {2p + ^ 

and where the coefficients and dim again depend on n,p,a, and 70 . 

Finally, one must ensure that the remaining equations coming from and the angular 
directions (%o = 0) are satisfied. This will require the following relation 


(A.27) 


hi) - 


= — 2Cr 


(A.28) 


This completes the analysis of the scalar sector. The remaining undetermined integration con¬ 
stants are thus: ,Crv,ci^\c^'^ ■ Note that the above functions reproduce the neutral case 

as 7 o — 7 > 0 . 
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Vectors of SO(p) 

The vector sector consists of 3p independent equations which correspond to the vanishing of the 
components: £ri,£vi and Mi. 


Constraint eqnations: The constraint equations are given by the Einstein equations = 0 
and are solved by 

diro = ro{l + N'yo)dyUi , (A.29) 

which are equivalent to conservation of stress-momentum. These are part of the conservation 
equations given by (1.6) in the rest frame. Similar to above we now proceed solving for the 
first-order corrections to the metric and gauge field under the assumption that the fluid prohle 
satisfies the above constraint (A.29). 


Dynamical equations: The remaining equations consist of p pairs of one Einstein equation 
= 0 and one Maxwell equation Mi = 0. The structure of these equations is the same as in 
the scalar sector. The Einstein equation £yi = 0 is schematically of the form, 


+ L^^’^\r)fyi{r) + LS"’^^(r)a'(r) = 5^i(r) , 
while the Maxwell equation Mi = 0 is, 

Ml^'^\r)a'l{r) + M^'"’^\r)ai{r) -h M{"’^^(r)/h(r) = Si{r) . 


(A.30) 


(A.31) 


We omit the expressions for the functions and M^, k = 1,..., 3. 

To decouple the system we differentiate £m once and eliminate all ai{r) terms in Mi. Doing 
so, one obtains a 3rd order ODE for /«i(r’) which can be written on the form 


dr 


r^+^/Cr) 


1 - Cl- 


dr 


l-Cl^l) 


fM 


= Syi{r) , 


Cl = 


V- 1 
1 -I- Njo 


70 


(A.32) 


It is possible to perform the first two integrations analytically and ensure regularity at the horizon. 
The first integration is straightforward while the second involves several non-trivial functions. The 
large r behavior of the fyi function is found to be 


fviir) 


.(1) 


- 1 - 


fjr) 

h{r)^ 


.(2) 


^nk 


- {dyUi)r + 


(^) I /q(^) 

(Kir + ^yi’rQ 


(A.33) 


k=l 


where the first two terms constitute the homogeneous solution and we find, in particular, that in 
order to ensure horizon regularity one must have 


q {2) ^ / 270(1 -b 70)( 0 ^-»^) -F ((9»7o) \ 

(l + 7o)^-^(l + iV7o) ) 


(A.34) 
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The remaining set of coefficients and are in general complicated expressions depending 
on the parameters in the problem. We therefore omit them as they provide no insight. Also, we 
notice that the sum in the function (A.33) vanishes in the neutral limit. 

Once the solution of /„ is given we can use £yi to determine Oj, 


aiir) 


.(1) 


+ 


.(2) 


h{r) 


+ E 

k=l 


y,nk 


-nk 


(k) 

a) r 


a(^) 

W ro 


(A.35) 


where the first two terms correspond to the homogeneous solution. Again, we choose to omit the 

coefficients and The remaining undetermined integration constants are thus: 

and 


Tensors of SO(p) 


There are no constraint equations in the tensor sector which consists of p{p + l)/2 — 1 dynamical 
equations given by 

£ij - = 0 . (A.36) 

This gives an equation for each component of the traceless symmetric perturbation functions fij^ 

A _ f \ M 

_ [r”+V(r)/b(r)] = -aijr^ l^2{n + 1) + h{r)^-^ , (A.37) 

with the same structure as Eq. A. 11. The spatial part of the shear tensor is 


The solution is given by, 




-Sijdku’^ 

p 


(A.38) 


fij{r) = Cij - 2aij (r^ - ^ (1 + 70) 2 log /(r)^ , (A.39) 

where horizon regularity has been imposed and the integration constant(s) Cij is symmetric and 
traceless. 


A.3 Solving the fundametally charged system 

As already explained, many of the differential equations appearing in fundamentally charged 
system are similar to the ones appearing for the Maxwell system. Instead of repeating these, in 
this section we provide the most important differences. 

Scalars of SO(p) 

The scalar sector consists of eight independent equations which correspond to the vanishing of 
the components: £yy,£rv,£rr,"£^£ij,£nn,£{ 4 >),Mai...ap+i and Mrh...ip with i G {n*}. 
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Constraint eqnations: There are two constraint equations; 8'^^ = 0 and M.rix...ip = 0 with 
i € {u*}. The two equations are solved consistently by 


dvro = -- 


ro(l + 270 ) 


-di: 


dvlQ = -- 


2n7o(l + 70) 


-diu^ . (A.40) 


n + 1 + (2 - n(iV - 2 )) 7 o ‘ n + 1 + (2 - n(iV _ 2 )) 7 o 

The first equation corresponds to conservation of energy while the second equation can be inter¬ 


preted as the charge density being constant in time. 


Dynamical equations: After the constraint Eqs. (A.40) have been imposed one is left with a 
system very similar to the one obtained in the presence of Maxwell charge {q = 0). It consists of 
the six equations given by the components: Syy, Syr, T^rSij, Snn, £’( 0 ), and M.ai...ap+i- 

The particular combination of TrSij and £( 0 ) gives Eq. (A. 11), but now with the relation 

r(r) = Tr/,,(r)- /^(r) , (A.41) 

(n -|-p -|- l)a 

and is solved by the same expression given by Eq. (A. 13). The equation for the trace Tr/jj is 
again similar to the q = 0 case and the solution can be put on the form given by (A. 15) with 
(A. 17), but where 


G{r) = - 


N{p -h 1) Tq 


n -|- p -|- 1 r 

and the coefficients are given by the expressions 
2pn‘^ f 2 (n + 1 ) + Cqo 


2 + 2 - 


Nn{p + 1) \ r\ 


n+p+1Jf 


r7o 


a = 


p + l V (n + 1)2 + C 7 o( 2 (n + 1 ) + Cyo) 


and j3 = 


P 


n 


p+l \n + 1 + Cjo 


(A.42) 


(A.43) 


The solution of Tr/^ dictates the perturbation of the dilaton held through Eq. (A.41). With 
Tr/jj determined, one can hnd the remaining perturbation functions as follows: fry from Syr and 


fyy from TrSij by a single integration while the gauge held perturbation a. 




can be obtained 


from A4 


vai...ap 


by a double integration. Einally, we note that remaining undetermined integration 


constants are equivalent to the = 0 case. 


Vectors of SO(p) 

The vector sector consists of 2>p independent equations which correspond to the vanishing of the 
components: Syi^ Syi and M.yrj^,,,jp_^ with j 7 ^ i. 


Constraint equations: The constraint equations are given by the Einstein equations £7 = 0 
and A4„r+...+_i = 0 with j / i. Eor each spatial index i one has a pair of equations that are 
solved by 


o ^ _ ^0(1 + 270) a 


2n7o(l + 70) o 

- j-ITt — IW—OvUi 


(A.44) 


l-(nV- 2 ) 7 o ‘ l-(niV- 2 ) 7 o 

The hrst equation corresponds to conservation of stress-momentum while the second equation 


censures that the charge density does not have any spatial gradients over the world-volume. We 


see that the current is more constrained compared to the case with Maxwell charge which is tied 
to the fact that the p-brane charge is not able to redistribute itself. 
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Dynamical Equations: After the constraint Eqs. (A.44) have been imposed the remaining p 
equations consist of 2 nd order differential equations = 0 of the form 

^ = Sri(r) . (A.45) 

Each equation can be integrated analytically in order to obtain the perturbation functions fm. 
Horizon regularity is ensured due to the form of the differential operator. We note that the 
homogeneous solution gives rise to two integration constants: and . 

Tensors of SO(p) 

It turns out that with the parametrization given by Eq. (A. 6 ), the equations for the tensor 
perturbations take the exact same form as found for the Maxwell charge given by the form 
(A.37). The solution is therefore, 

fij{r) = Cij - 2aij (^r* - ^ (1 + 70 )^ log f{r)j , (A.46) 

with given by Eq. (A.38) and regularity of the horizon has been imposed. The constant Cij 
is again symmetric and traceless. Because of this closed-form expression, the shear viscosity will 
take the same form as found for the system with Maxwell charge. 


A.4 Fixing the integration constants 

We have determined the first-order derivative corrected solution for both types of branes. Both 
solutions are completely determined up to a set of integration constants, , (^\ ci‘\ 

Crv, Cij (and for g = 0). As explained in Sec. 1.1, these constants are hxed by virtue of 
asymptotic flatness and choice of fluid frame (gauge). The latter freedom exactly corresponds to 
constant 0{d) shifts of the parameters of the 0{d^) helds, i.e., it parameterizes the homogeneous 
solution to the above differential equations. Indeed, for the zeroth order Maxwell solution (A.4), 
consider order 0{d) constant shifts tq — >■ tq -|- (5ro, 70 —>■ 7o + <^7o and a constant gauge shift 
—)■ + Stty. Moreover, by redefining the r coordinate as. 


rr ( 1 -h 7o(n(51ogro(51og7o)G(r)) , 


(A.47) 


with G(r) given by (A.16), the angular directions do not receive hrst-order contributions in accor¬ 
dance with the gauge choice (gd)nQ = 0. The resulting expressions after the shifts and coordinate 
transformation are of course still solutions and the overall change exactly corresponds to the 
homogeneous solution to the above differential equations in the scalar sector. More precisely, one 
can relate the integration constants to the shifts by, 


= -2p(nSlogro dlogyo) , 4 ^^ = 


H) - 


Sd'}} 


ro\/Njo(l +jo)' 
70 


= — ndlogro- ^ (51og7o- 7 _ 

2(1 + 70) 4^70(1 + 7 o) 




(A.48) 
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For the vector sector one finds that the homogeneous part of the above solution corresponds to 
global constant shift of the boost Ui ^ Ui + 5ui and in the gauge Oj —>■ + dai. With the same 
choice of radial coordinate, one has 




db _ 


6(Xi 


c- = — 




(A.49) 


Similar expressions relating the shifts to the integration constants can be derived for the fun¬ 
damentally charged solution. In the following we require all the 0{d) shifts to vanish. In the 
effective fluid description this exactly corresponds to choosing the Landau frame. This fixes the 
integration constants, = 0 (and = 0 for g = 0). 

To fix the remaining integration constants, we now impose asymptotic flatness to 0[d). In 
order to do so, we must first transform our results back into Schwarzschild form. In addition 
to making the asymptotics more transparent, this is also needed for extracting the effective 
hydrodynamic currents. In order to change coordinates, we use the inverse of the transformation 
stated in equation (A.l) to 0{d). The inverse transformation is worked out iteratively order by 
order. To first order, the transformation from EF to Schwarzschild coordinates is found to be. 


V = t + r^ + 


{t r*) {droT^dtro + d^^r^dao) -h x* {droV^^dirQ 


+ 0{d^) , 


a=x' + 


{t r^)dtu!‘ + a^dju" 


r^ + Oid^) . 


(A.50) 


It is now possible to express all the helds in Schwarzschild coordinates and impose asymptotic 
flatness. This leads to 

Crv = 0 , = 0 , = 0 , Cij = 0 , (A.51) 

for both types of branes. Having obtained the full first-order derivative corrected asymptotically 
flat solutions for both the Maxwell and fundamentally charged brane, it is now possible to read 
off the effective hydrodynamic currents using standard methods (see also [9]). Our results for the 
transport coefficients are presented in Sec. 1.2 and 2.1. 
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